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A simple phenomenological model of a binary granular mixture is developed and investigated 
numerically. We attempt to model the experimental system of where a horizontally vibrated 

binary monolayer was found to exhibit a transition from a mixed to a segregated state as the filling 
fraction of the mixture was increased. This model is found to reproduce much of the experimentally 
observed behaviour, most importantly the transition from the mixed to the segregated state. We use 
the model to investigate granular segregation mechanisms and explain the experimentally observed 
behaviour. 
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I. INTRODUCTION 

Granular systems exhibit a wide range of intriguing 
and often counterintuitive phenomena. Segregation of 
two or more species of grains due to vibration or shear- 
ing is one such example |3, 3- Many mechanisms, includ- 
ing buoyancy, temperature gradients, differing angles of 
repose, and differing roughness of particles, have been 
proposed 0, 0, IE • Recently, segregation of a vibrated 
binary mono-layer has been demonstrated in a series of 
experiments that imply/suggest the existence of 

a segregation critical point as the filling fraction is varied, 
with associated growth of fluctuations and timescales in 
the vicinity of that point. 

In this paper we propose a phenomenological model 
of the experimental system which captures its essential 
features. We show that there is qualitative agreement be- 
tween the model and the experiment and, inparticular, 
that the quantitative measures reported in [ij, |3| are of 
the same form in both cases. We then use the model to 
study segregation mechanisms, giving one definite mech- 
anism and demonstrating that a second may also play a 
role. 

The experimental system used in 0, Q consists of a 
smooth horizontal tray of dimensions = 180mm by 
Ly = 90mm vibrated sinusoidally parallel to its major 
axis at a frequency of / = 12Hz and amplitude oi A ~ 
1.74 ± 0.01mm. The grains are high-precision phosphor- 
bronze spheres (denoted by c in the following) of radius 
Rc = 0.75mm and mass mc — 16.8fig and poppy seeds 
(denoted by p in the following) of average radius and mass 
Rp = 0.54mm and m^ = 0.52/ig respectively. The poppy 
seeds are rough, non-spherical polydisperse [1| particles. 

Particles placed on the oscillating tray predominantly 
move periodically with the driving. They also have some 
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quasi-random component to their motion, caused in the 
case of the ps by their non-sphericity which means that 
they do not roll smoothly but instead 'slick and slip' and 
scatter. The cs are affected less by the driving since they 
roll as the tray moves beneath them. However they too, 
when placed individually on the tray, are seen to have 
some quasi-random component to their motion. The 
overall impression is of a granular bed of agitated par- 
ticles with the predominant motion being parallel to the 
driving. 

The control parameter used in the experiments was 
the 2-dimensional filling fraction, termed the compacity 
C = {NcirRl+NpirRl) / {L^Ly), with A^^ fixed. The value 
of Nc is also important but for simplicity the experiments 
focussed on Nc ~ 1600 and varied Np. The system was 
quasi 2-dimensional in that 0.49 < C < 1.12, high C be- 
ing achieved by poppy seeds 'riding up' and overlapping 
each other to some extent, although never so much that 
they were above and overlapping the cs. 

Starting from a uniformly mixed initial distribution, 
the final state reached by the system varies with C. 
For C < 0.65 the system remains mixed whilst above 
that value small mobile segregated clusters of cs form. 
As C increases, these clusters grow in size and become 
anisotropic, forming stripes perpendicular to the direc- 
tion of driving for C large. For C > 0.93 the cs within 
these stripes crystallise to form a densely packed hexag- 
onal lattice. These three phases were termed binary gas, 
segregation liquid, and segregation crystal. The existence 
of a phase transition from the binary gas to segregation 
liquid, with an associated critical value of C, was re- 
g^ed in 0. The details of these results are given in 

The experimental system has several appealing prop- 
erties from the point of view of studying granular mat- 
ter: except at high compacities where stripes form, the 
final state reached is a function only of the compacity, 
i.e. the initial conditions are not relevant. The particles 
are always in contact with the tray and hence are al- 
ways effectively 'thermalised'. This is in contrast to the 
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behaviour of particles in granular systems such as sand- 
piles where the grains spend most of their time locked 
in position, or some vertically vibrated systems where 
much of the time is spent in free flight The constant 
'thermalisation' of the particles, their ability to explore 
many possible states, and the irrelevance of initial con- 
ditions suggests that the system might be a good choice 
for studying the 'statistical mechanics' of granular mat- 
ter. However detailed balance, equipartition of energy 
and other rigourous features of equilibrium statistical me- 
chanics are of course not obeyed in this system [To|. 

II. DESCRIPTION OF MODEL 

We here develop a phenomenological model of the ex- 
perimental system. The aim is to capture the essential 
features of the experiment in the model, thereby discov- 
ering what those features are and in particular what the 
segregation mechanism(s) is(are). 

The following features are, we believe, necessary: there 
is a tray of dimensions x Ly whose base moves with 
sinusoidal velocity Asva(ujt). There are two particle 
species, named c and p. The particles moving relative 
to the tray surface are subject to a frictional force. The 
particles feel a randomisation of their velocity caused by 
their stick-slip motion with the oscillating tray and their 
non-sphericity (thermalisation). The particles collide in- 
elastically. 

We therefore model the system as 2-dimensional with 
the particles behaving as hard disks of mass ttIq and ra- 
dius i?Q,, where a denotes the species (c or p). Except 
during collisions the particles obey the Langevin equation 

ma-Vai = -lai^ai ~ ^tray) + Vatit) (1) 

where vtray = sm{Lut) and 7 provides a linear damp- 
ing, tjalt) is Gaussian white noise of mean zero and 
standard deviation {rjait) ■ rjait')) = 2a^6{t — t') and 
this provides the 'thermalisation'. The particles interact 
through smooth hard-disk inelastic collisions with coeffi- 
cient of restitution r^^f}, i.e. in the centre of mass frame, 
V||i — r^jVp, v±i vxi, where V|| and vj_ are the ve- 
locity components parallel and perpendicular to the line 
joining the centres of the particles i and j. For simplic- 
ity, the disks have been taken to be smooth sided so that 
angular momentum can be ignored. A similar model has 
been used in ^ to describe colloidal particles driven by 
an external electrical field. Also IT^ have used a similar 
model without noise to model granular particles driven 
by a vertically oscillating air column |T^ . 

These are the essentials. We have also kept the walls of 
the box stationary for simplicity but modelled the motion 
of the end walls by considering that in collisions they 
have a velocity max(iAi sin(wt), 0) for the left wall and 
mm{iAi sm{uit) , 0) for the right wall. We felt that this 
was necessary in order to model the low-density region 
near each end wall caused by the vigorous collisions with 



the end walls. The width of the low-density region is 
independent of the system size, thus we used Ai = A/10 
rather than Ai — Ain the results described here in order 
to reduce the 'finite size effect' of this region as the system 
size is changed. 

There are many approximations of the real system 
made here, the most significant ones are: the friction 
term 7V is only an approximation, it is chosen as be- 
ing the simplest possible form. The noise is in fact due 
to the stick and slip interactions between the particles 
and the oscillatory surface, and also their non-sphericity 
when interacting with each other and with the tray. We 
do not try to directly model this since it would require 
detailed specification of the shape of each particle and 
its actual interaction with the tray, which is not known. 
Even a single high-precision phosphor-bronze sphere con- 
ducts a quasi-random walk when placed on the oscillating 
tray, indicating that the randomness can depend on very 
small imperfections of the particles and the tray (and 
possibly also in the driving). Both because of this im- 
mense difficulty and in order to have a reasonably simple 
model whose behaviour we can understand, we instead 
choose to include the noise phenomenologically. We as- 
sume that the noise the particles receive is independent 
of their neighbours and of the phase of the tray cycle, 
both of which are unlikely to be accurately what hap- 
pens in the experiments. The assumption that the noise 
is Gaussian and white is an approximation. For simplic- 
ity we are using a 2-dimensional model, which ignores 
the overlapping of ps and the rolling of particles. The 
particles and walls are assumed to be smooth and so an- 
gular momentum is not considered. We have ignored the 
polydispcrsity of the ps. Polydispersity was included to 
check its importance and was found to leave the qual- 
itative behaviour unchanged. The final approximation 
is that the coefficients of restitution are constant which 
is a commonly made one |l4j . Despite these simplifica- 
tions and approximations, in Section llTTI we show that our 
phenomenological model captures much of the behaviour 
observed experimentally. 

A. Parameter Values 

Static parameter values such as mass and size can be 
measured reasonably accurately, for the poppy seeds we 
have used the mean values of a sample of measurements 
|l5j| . Dynamic parameters were less accurately known, 
7q, was estimated from the distances travelled by single 
particles striking the moving end walls. Using the result 
for the noiseless case, 

x{oo) = x(0) +v^{0)m/-f (2) 

and estimating VxiO) to be equal to the maximum veloc- 
ity of the end wall gives approximate values for 7^ . The 
value for the noise is the hardest to determine since no 
velocity or accurate r{t) path measurements were avail- 
able. We merely estimated that the mean square velocity 
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due to the noise, (v^^) = a'^ma/^fa should be equal to 
~ (A/F)'^ where the factor F = 3 ior the p case and 
F = 13 for the c case. 

Clearly these last estimates are rather crude. However, 
extensive study of a wide range of these parameters has 
shown that the qualitative features are robust to varia- 
tion of these estimates. 

The coefhcients of restitution, ra^fs — Tp^a. are esti- 
mated to be: rc,c = 0.9, rc,p = 0.2, Vp^p = 0.1, Tc.w = 0.9, 
'i'p,w = 0.2 where w denotes a side wall. Other values 
have been studied, but the qualitative behaviour is un- 
changed. The main effect of increasing(decreasing) r^^^ is 
to increase(decrease) the granular temperature which in 
general merely moves the onset of segregation to slightly 
higher (lower) compacity values. 

The parameters then are the following. 



property 


c value 


p value 


nia mass 

Ra radius 

7q, damping term 

Uxa noise term 

Uy noise term 

ra.,/3 


1.6800 X 10"^ 
7.5000 X 10-"^ 
4.3636 X 10-*^ 
8.1819 X 10"^ 
8.1819 X 10-** 

^c,c — '^c^w — 0.9 

rc,p - 0.2 


5.2000 X lO"'^ 
5.4000 X lO-"* 
1.0000 X 10^5 
1.0000 X lO"'^ 
1.0000 X lO"'^ 
^c,p — '^p^w — 0.2 
rp.p = 0.1 



Table III Al The parameters used for the simulation re- 
sults described below (SI units are used at all times unless 
stated otherwise), ax a and are the x and y compo- 
nents of a. 



B. Simulation Method 

We have simulated the model via an event driven code 
[T^ . The process is as follows: 

1. For each particle, predict when it will next collide. 

2. Identify the first collision to occur. 

3. Move the particle(s) involved so that they touch. 

4. Update the velocities of the particle(s) involved 
(change in velocity due to damping and noise since 
they were last updated). 

5. Colhde the particle(s) p^ . 

6. Re-predict the next collision(s) of the particle(s) 
and their neighbours. 

7. Repeat from 2 until time has advanced by 

^minupdate / 

8. Update all particles that have not been updated in 

the last tminupdatef^ SCCOuds. 

9. Repeat from 1 imtil time has advanced by Uakedata- 

10. Record data. 

11. Repeat from 1. 

The prediction of collisions assumes that the particle's 
velocities do not change during their motion. The er- 
ror caused by this is small provided that tmmupdate ^ 



min(TQ) where Tq = rUa/^a is the time constant of the 
velocity decay. We set tnunupdate = 0.01 x min(Ta). 

The hard sphere model with inelastic collisions can un- 
dergo inelastic collapse "vf] , where the particles undergo 
an infinite number of collisions in a finite time. Clearly 
a simulation that implements each collision will 'stall' 
in these circumstances. One way around this unphysi- 
cal singularity is the tc model jig which prevents the 
collapse by setting Va^p — 1 for any particle which col- 
lided within the last ttcoimin seconds. We found that for 
ttcoimin < 1 X 10"'* the results did not change. For the 
results presented here we chose ttcoimin — 5 x 10~^. 

The initial conditions were created by running a re- 
duced size system Nc/25, iVp/25, of twice the aspect ra- 
tio with particles initially placed randomly on a square 
lattice. An external force was applied to compress the 
particles into an area of size L^/Si x Ly/5. The particles 
were then allowed to move in this box without the ex- 
ternal force until equilibrated. The initial condition was 
then created by tessellating the full system size with 25 
replicas of the reduced system. To prevent long-range or- 
der the replicas were randomly inverted in both the x and 
y directions (the differing noises received by the particles 
in conjunction with the chaotic behaviour of the particles 
would rapidly remove any correlations in any case). Dur- 
ing this stage, both particle species had the same prop- 
erties, including Va^p = 1, except for their radii. This 
method produced initial conditions which appeared to 
be as homogeneous as those for the experiment. 

The simulations were run on standard PCs. For com- 
parison with the experiments our results are for the same 
system size (18cm x 9cm), except where stated otherwise. 
The aspect ratio is here 2 : 1 in all cases. 



III. RESULTS 

As with the experiments, we used the number of ps, Np 
as our control parameter. The main quantities measured 
were those found experimentally 0, S 01 and related to 
the cs only: the mean stripe width (MSW), and the local 
density, pi, of the ith particle. We also visualised the 
system and watched its behaviour. In addition, we also 
measured the area available to the ps, and the kinetic 
energy or 'granular temperature' of the particles. 

The MSW measures the mean width of the stripes in 
the X direction by pixellating an image of the system, 
deciding which pixels Px^y are within a c domain, then 
running along each pixel row Py and counting the width 
and number of domains. All the rows Py are summed and 
the average domain width found. Whether Px.y is within 
a domain is decided by blurring the image with a gaussian 
smoothing function and then setting a threshold. This 
measure was used in [ij, even when the domains had 
not formed into stripes since it provided a simple measure 
of the domain sizes in the longitudinal direction. 

The normalised local density pi was found by Voronoi 
tessellation of the cs, such that each Ci has around 
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FIG. 1: The time evolution of a simulation showing the coars- 
ening into domains for a system of 1 /4 the area of the exper- 
iment and C — 0.6721. The times are, top left: t=0.04s, top 
right: t=4.18s, middle left: t=8.37s, middle right: t=16.75s, 
bottom left: t=33.51s, bottom right: t=62.83s. The ps are 
coloured black. 



it a polygonal area all points of which are closer to Ci 
than to any other c. pi is then the minimum possible 
area, which is 2-\/3i?ci divided by the polygonal area. 
Polygonal cells on the edge of the system which are not 
bounded are discarded. 

In addition to the MSW as a measure of the amount 
of coarsening, we also measured the area available to the 
ps. This 'available area' is just the fractional area of the 
system in which a p could be placed without overlapping 
a c. Each c has a circular 'excluded area' around it of 
radius Rc + Rp, inside which the centre of a p cannot be 
placed. If all Nc cs are widely separated, the available 
area is 1 — Nc tt{Rc + Rp)^ / LxLy, whilst if all the cs 
are hexagonally packed in one domain the available area 
will be larger (w 1 — Nc 2\/iRc^ / LxLy) since the ex- 
cluded areas now overlap. Thus the available area gives 
a measure of how segregated the cs are. In Section IIV Bl 
we will show that the available area is relevant to noise 
segregation. 

Figure ^ shows 6 images of the evolution of a coars- 
ening system. It can be seen that the initial segregation 
into relatively small domains is rapid. This is followed 
by slower coarsening as domains merge and as larger do- 
mains grow at the expense of smaller ones that 'evapo- 
rate'. This growth rate is also clearly seen in the time- 
series plots of figures 13 and 21 Figure [21 shows three im- 
ages of the system at differing compacities for late times, 
the lowest compacity (C=0.446) shows a binary mixture, 
the second (C=0.582) a segregated 'liquid' which has mo- 
bile and transient clusters and shows a slight anisotropy, 
and the third (C=0.717) shows a system that has coars- 
ened into stripes. The first two are in a steady state 
whilst the third is still evolving slowly due to the high 




FIG. 2: Examples of a binary mixture, segregated liquid, and 
segregated stripes. Compacities are 0.446, 0.582, and 0.717, 
respectively, ps are coloured black. The pictures were taken 
after 100 seconds for the gas and 200 seconds for the liquid 
and crystal, by which time the first two have reached a steady 
state whilst the third is still slowly evolving. 



compacity (for a strictly 2-dimensional system) which 
causes particles to be 'blocked' by other particles. 

Figure shows MSWc(t) for several compacities. As 
C increases, the cs, which were initially mixed with the 
ps, coarsen into domains whose size increases with C. 
The early-time coarsening is rapid, followed by slower 
coarsening and then saturation at some relatively steady 
value. For the lowest compacities, there is no coarsening, 
the system remains in a mixed, disordered state. The 
highest C value shown has not reached a final steady 
state, the stripes are still moving and merging at a very 
slow rate compared to the initial coarsening. At higher 
compacities, the system becomes blocked or jammed, 
the particles being unable to rearrange themselves in 2- 
dimensions, the slow timescales of the top curve show 
the onset of this jamming. This jamming does not oc- 
cur to the same extent in the experimental system since 
the particles do not always form a monolayer and ps can 
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FIG. 3: Plots of MSW against time for various compacities. 
From the right, top to bottom, the compacities decrease from 
C=0.717 to (7=0.401 in uniform increments. Thus it can be 
seen that for low C the system does not coarsen whilst when 
C is increased, the system coarsens to a roughly constant 
value which increases with C. For the largest value of C, 
the system has undergone rapid initial coarsening but then 
slowed as the large domains move more slowly, especially at 
this high compacity where the system is becoming somewhat 
'jammed'. The MSW was measured every 0.0209 seconds. 



move out of the way of cs by 'riding up' on top of each 
other. 

Figure El shows a plot of the value of the area available 
to the ps against time, it is similar to that for the MSW. 
It is a more repeatable measure than the MSW since 
it is less affected by stripes merging, for example, two 
runs with the same parameters naturally differ due to the 
chaotic behaviour of the particles. This means that for 
high compacities merging of stripes in two realizations 
at the same parameters may occur at different times, 
causing the MSWc{t) curves to differ between runs at 
late times (since MSWcit) is inversely proportional to 
the number of domains). The available volume is less 
affected by this and thus provides a 'cleaner' measure of 
the coarsening. Thus for later results we will use the 
available volume, although the MSW provides similar, if 
more noisy, results. Experimentally, the available volume 
has not been measured since the experimental system 
was only imaged in it's central region, thus the number 
of cs changes with time making the measure somewhat 
arbitrary (the changes due to cs entering or leaving the 
system would be more significant than the coarsening). 
We note that for the highest compacity curve, the area 
available has increased by « 0.15 whilst the maximum 
range of the available area stated before is « 0.30. 

In figures and [HI we plot the late time values of the 
MSW and available volume as a function of the compac- 
ity. This was done by fitting exponentials to the time- 
series of the type shown in figures |21 and 01 and using 
the late-time values. We do not claim that the time se- 



FIG. 4: Plots of the available volume for the ps against time 
for various compacities. From the right, top to bottom, the 
compacities decrease from C=0.717 to C=0.401 in uniform 
increments. The conclusions are similar to those for figure |3 
except that the curves here always remain in the same order. 
The available volume was measured every 0.838 seconds. 



ries are exponential, but the fits provide us with a rea- 
sonable measure of the late-time values. The MSW vs 
C curve is the same measure as that used experimen- 
tally in (3, ^2d\ to claim a mixed state to segregated state 
phase transition for the experimental system. In a 
square root curve was fitted to the right hand side of the 
data, the data to the left of the transition being taken 
to be roughly constant, i.e. MSWsat{C) — B for C < 

CtransMon and MSW,at{C) = B + D^C^t ransition 

for C > Ctransitiom wherc B and D are constants. From 
figure Owe see that this is not the case for the simula- 
tion data, which, as C is increased, initially rises slowly, 
then increasingly rapidly before slowing again for large C. 
There is no indication of a discontinuity in the gradient 
of the order parameter, merely a rapid increase. In "^o"! 
the same experimental data is presented with a sigmoid 
shaped 'guide to the eye' curve rather than a square root. 
In our opinion the simulation data is similar to the exper- 
imental data but not to the square root form suggested 
in ,2i|. The main conclusion from both the experimental 
and simulation MSWsatiC) curves is that MSWsatiC) 
has a roughly sigmoid shape, with MSWsat{C) increas- 
ing rapidly with C in the central region. 

The simulation value of C at which the rapid increase 
occurs (« 0.58) differs somewhat from the experimental 
value (0.647 ± 0.049) 0. In the next section we show 
how this value changes continuously as we vary the noise 
strength or other parameters. Thus, by increasing the 
noise strength we can increase the value of C at which 
the rapid increase occurs. 

Figure [7| shows histograms of the local Voronoi density 
for various compacities, the results are qualitatively simi- 
lar to the experimentally reported ones 1] . At the lowest 
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FIG. 5: Plot of the saturated (late time) mean stripe width 
as a function of compacity. The saturated values were found 
by fitting exponentials to the MSW(t) curves. Two runs were 
done at each compacity, and the results averaged. The error 
bars are based on the difference between these two runs. 
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FIG. 6; Plot of the saturated (late time) area available to 
the ps as a function of compacity. The saturated values were 
found by fitting exponentials to the available-area(t) curves. 
Two runs were done at each compacity, and the results av- 
eraged. The error bars are based on the difference between 
these two runs. Note that the error bars are smaller than for 
the MSW vs C plot since, as stated in the text, the available 
area is a 'cleaner' measure. 



C values, corresponding to the mixed state, the distri- 
bution is peaked at low densities as one would expect 
for unclustered cs. At high C the distribution is peaked 
at large densities as one would expect for clustered (i.e. 
segregated) cs. There is a crossover between these two 
cases, with figure E|i showing a broad histogram due to 
almost the full range of densities being present in almost 




FIG. 7: Normalised histograms of the Voronoi densities plot- 
ted against normalised density, for C=0. 400728(a), 0.559(b), 
0.582(c), 0.604(d), 0.627(e), and 0.740(f). Note the crossover 
with increasing C from a single peak on the left to two peaks, 
then a larger peak on the right. For each compacity the data 
was measured in the steady state from 45 frames of 1600 cs 
each. 



equal weights. Unlike the experimental results however, 
the high C distribution also has a peak at low density 
caused by a small fraction of isolated cs which are not 
present in the experiment. There is a crossover between 
these two extremes as C is increased, the central 4 fig- 
ures (b-e) show this crossover. Following 0, we plot the 
location of the peak(s) and their widths as a function of 
compacity in figure|Sl In I] the peak width used was the 
full width 3/4 maximum since the peak did not extend far 
enough above the rest of the distribution for a full width 
1/2 maximum to be meaningful. Here we did the same 
although a full width 3 /4 maximum also does not exist in 
one case. Figure |H1 is qualitatively different from its ex- 
perimental equivalent and we conclude that although the 
histograms are qualitatively similar, the results derived 
from them are not. 

The results presented here show that our model repro- 
duces the segregation and its qualitative behaviour with 
C of binary gas, segregated liquid, and segregated stripes 
as seen experimentally. This behaviour is not immedi- 
ately obvious from the model rules - it emerged from the 
set of rules which we believed contained the important 
microscopic features of the experiment. We have also 
shown that the MSW as a function of time and its satu- 
rated value as a function of compacity behave in a quali- 
tatively similar manner to the experiment. These repro- 
ductions of experimentally observed behaviour lead us to 
conclude that the 'necessary features' listed in section UTI 
capture the essential behaviour of the system. Nonethe- 
less, our phenomenological model does not quantitatively 
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FIG. 8: Plots of data derived from the Voronoi histograms as 
shown in figure[7| Squares show the position of the peak(s) in 
the histograms, representing the 'most probable' densities, the 
lower curve is for the leftmost peak and the upper curve is for 
the rightmost peak. Filled circles show the full width of each 
peak at 3/4 maximum (where the peak does have a distinct 
3/4 maximum). The curves differ from their experimental 
equivalents where there was only one maximum (thus just two 
single curves) which moved steadily to the right and where 
the width was peaked at an intermediate compacity. These 
differences are due mainly to the presence in the simulation 
density histograms of a second peak at low density which is 
always present and is caused by a small number of isolated 
cs, something not observed in the experiment. 



reproduce the experiment (it was never expected to do 
so), in particular the Voronoi histograms differ from the 
experimental ones sufficiently that the data derived from 
them (figure IHll is qualitatively different. 



IV. SEGREGATION MECHANISMS 

Having shown that our model is relevant to the exper- 
iment, we now use it to investigate segregation mecha- 
nisms. 



A. Segregation due to Oscillatory Driving 

That the domains are anisotropic and indeed form 
stripes at high C indicates that the anisotropy of the 
driving is significant. Experimentally the stripes form 
perpendicular to the driving even for aspect ratios greater 
than 1, e.g. Ly/L^ = 2. 

The side-to-side driving causes the two species to move 
at different rates due to their differing masses and fric- 
tion coefficients. A single particle will oscillate with 



{xa{t)) = A/ -y/f -I- uP-(mal^aY sin(u;t) where the aver- 
age is over the noise. Thus {xc{t)) = 6.00 x f sin(wi) 
and {xp{t)) — 4.30 x f 0^* sin(wi). Thus the ps would 
'like' to move a distance of order their radius during a 
cycle whilst the cs hardly move. Consider a state of only 
ps, if we remove the noise then all the ps move in the 
same sinusoidal way, if we transfer to the (non-inertial) 
reference frame in which they are at rest, we see that this 
is identical to the state with no sinusoidal driving (apart 
from edge effects at the walls) for which the dissipation 
causes the particles to be stationary. The same can be 
said of a state of only cs. In a mixed state, however, the 
ps will collide with the cs and the system will 'scatter' 
into a different state. Stable states, i.e. those that do not 
undergo further scatterings, will be those for which the cs 
are separated from the ps in the x direction by distances 
of at least the amplitude of the ps' oscillations. For the 
packing fractions considered here, this can only happen 
by the two species segregating into domains. The fact 
that the area available to the ps increases with time is 
consistent with this interpretation. This argument only 
holds when the driving is not too large compared to the 
dissipation. For example if the amplitude of oscillation 
were of order the system size then it would not be valid. 
An argument similar to this was also given in jl^ . 

We have studied systems with very low noise and the 
results obtained agree with the heuristic argument given 
above. The addition of noise, which causes particles to 
diffuse, will tend to cause mixing. Thus there is a com- 
petition between the periodic driving which causes seg- 
regation and the noise which prevents it. This is shown 
in the results presented below. 

For the standard parameters but with ctq O.I x CTq,, 
the system segregates for all the compacities studied in 
section Hill (see figure indicating that, as expected, 
the noise acts to prevent segregation. To confirm this 
we then gave the two species identical parameters (the p 
parameters) and set all r^.^ = 1. The only difference was 
that the ps experienced the periodic driving term jvtray 
whilst the cs did not. The results are shown in figure [TUI 
and demonstrate that a difference in the periodic driving 
alone can cause segregation. 

For C — 0.498 with the same parameters except 7q, — > 
2 X 7q (to reduce the granular 'temperature' which is 
higher than usual due to the r = I) we varied a and, as 
shown in figure ITTI the system segregates for low a and 
remains mixed for high a. 

These results confirm that the state of the system is 
a result of a competition between the periodic driving 
which causes segregation and the noise or 'granular tem- 
perature' which prevents it. Varying our control parame- 
ter, Np, causes us to go from a mixed to a segregated state 
because increasing Np both reduces the granular temper- 
ature (since it increases the number of collisions which are 
highly inelastic) and also increases the 'pressure' that the 
cs feel due to collisions with a greater number of oscillat- 
ing ps. Figure IT^ shows results for all properties set to p 
values including all r = O.I but with only the ps feeling 
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FIG. 9: Plot of the area available to the ps at late times as a 
function of compacity. The curve shows standard parameters 
but with low noise (ctq ~» 0.1 *(Tq) results, the two box points 
are the equivalent results for the standard parameters. Thus 
it can be seen that all compacities have segregated. For high 
C a few large domains form whilst at lower compacities there 
is enough space for a larger number of small domains to be 
stable. Note that the area of the system is 1/4 that of the 
experimental system. The late-time results were found by 
fitting an exponential to the available area vs time curves. 



FIG. 10: Plot of the area available to the ps at late times as 
a function of compacity. The data is for a system where both 
particle types are the same (standard p parameters) except 
that only the ps feel the periodic driving. We get a mixed 
state for low C and a segregated state for high C, as we did 
for the standard parameters. Note that the area of the system 
is 1/4 that of the experimental system. The late-time results 
were found by fitting an exponential to the available area vs 
time curves. 



the periodic driving, showing that the results remain of 
the same form for r < 1. 

Finally, we find that the system with standard parame- 
ters does not segregate if the oscillatory motion is turned 
off. 

These results show that the differential oscillatory driv- 
ing can cause segregation and present good evidence that 
it is the responsible mechanism in the system in conjunc- 
tion with the noise which acts to prevent segregation. 
The transition from a mixed to a segregated state as C 
is varied is due to the compacity changing the relative 
strengths of these two competing effects. 

B. Noise Segregation 

It was suggested in 0, that the segregation mecha- 
nism might be similar to the depletion interaction in equi- 
librium binary systems |2]| . As an example consider a 
colloidal suspension containing non-adsorbing polymers. 
In the ideal case where there are no forces present, the 
free energy depends only on the entropy of the system. 
Treating the polymers as spheres 22] , it is clear that each 
colloidal particle has an 'excluded volume' around it of 
radius Rcoiioid + Rsphere which the centre of the polymer 
cannot enter. Thus the volume available to the polymers 
is the volume of the system less the excluded volumes 
around the colloidal particles and the system edges. How- 
ever, the excluded volumes overlap when colloidal parti- 
cles are closer than 2{RcoUoid + Rsphere), thus the vol- 
ume available to the polymers, and hence their entropy. 
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FIG. 11: Plot of the area available to the ps at late times as 
a function of the noise strength a /a standard, for C — 0.498. 
Increasing the noise strength brings us from a segregated to 
a mixed state. Note that the area of the system is 1/4 that 
of the experimental system. The late-time results were found 
by fitting an exponential to the available area vs time curves. 



is larger if the colloidal particles are close to each other. 
This entropic 'effective potential' can be large enough to 
cause the colloid to coagulate. This mechanism was one 
reason for measuring the area available to the ps in our 
simulations. Segregation has been observed in simula- 
tions of hard spheres of two different sizes |2^ for large 
size ratios, e.g. R1/R2 = 10, and also experimentally in 
binary mixtures of hard sphere colloids [2J| . 
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FIG. 12: Plot of the area available to the ps as a function 
of compacity, the parameters are as for figure 1101 except that 
here the coefficients of restitution are all 0.1 rather than 1. 
It can be seen that this does not qualitatively change the 
results. Note that the area of the system is 1/4 that of the 
experimental system. 



We may equivalently view the entropy argument from 
the kinetic point of view. Two particles which are close to 
each other such that no third particle may fit in the space 
between them will feel a pressure on all sides due to col- 
lisions with other particles except on their neighbouring 
sides. Thus the particles feel an effective attractive force. 
This pressure argument may be extended to systems that 
are not in equilibrium, for example the granular experi- 
ment studied here. The size ratio Rd Rp is much closer 
to unity than for simulated equilibrium segregating sys- 
tems pJI, implying that the difference in size alone does 
not cause segregation. 

In our out-of-equilibrium system there are several pos- 
sible differences between the two species besides a dif- 
ference in size. For our system, it seemed possible that 
the lighter, faster moving ps might, through the differen- 
tial pressure mechanism, cause the cs to coagulate even 
in the absence of periodic driving. As stated before, 
this was not observed for the standard parameters. We 
therefore increased the noise of the ps and/or reduced 
the noise of the cs in order to increase the p to c tem- 
perature and hence pressure ratio. Noticeable segrega- 
tion occurred over a wide range of compacity values for 
CTj, — > 10 X (Tp, for (Tp — > 2 X (jp and CTc ^ 0.1 x a^, and 
also for CTc — > 0.05 x ac- Whilst the first of these cases is 
outside the reasonable parameter range, the second and 
third are at parameters which might be physical. Fig- 
ure El shows an example for the third case. Notice that 
many of the cs have coagulated at the walls as one would 
expect since the pressure argument for two-particle at- 
traction also applies to the particle- wall case. Although 
this marked congregation at the walls is observed in ex- 
periments with colloids, it is not observed in the exper- 
iments of P, 13 ■ Re-introducing the periodic driving of 




FIG. 13: Segregation of cs for the standard parameters but 
with no periodic driving and ctc — > 0.05 x tJc. Notice the 
domains of cs at the edges of the system. Note that the area 
of the system is 1/4 that of the experimental system, ps are 
coloured black. 



the end walls prevented congregation at the x — Q^L^ 
walls whose large momentum transfer to the particles, 
as stated earlier, gives rise to low-density regions next 
to them. However, for all parameter values studied that 
displayed segregation with no driving, the cs still congre- 
gated at the y — 0, Ly walls when the driving was turned 
on. Since stripes that touch the top and bottom walls are 
stable in the experiment, it seems unlikely that agitation 
due to the motion of the top and bottom walls is what 
prevents the liquid domains from coagulating there. It 
is possible to remove this experimentally unobserved ef- 
fect by giving ax and ay differing values such that axe 
is lower and aye is higher than that needed to produce 
segregation. Whilst we had previously kept ax = ay for 
simplicity, it is reasonable that ax > ay since random mo- 
tion caused by sticking and slipping is likely to be larger 
in the direction of driving. This extra modification pro- 
duces segregation without cs congregating on the walls 
(provided that the x — 0, Lx walls are 'driving') whether 
there is periodic driving of the tray base or not. Figure 
1141 shows an example. 

We therefore conclude that this differential pressure 
segregation mechanism may play a role in the experi- 
ment. We had to 'tune' the parameters which implies 
that the mechanism is less robust than the oscillatory 
driving mechanism. Accurate experimental measure- 
ments of the parameters would help to resolve whether 
this mechanism is indeed present. Directly distinguishing 
the two mechanisms mentioned would require accurate 
tracking of all particles and their collisions, coupled with 
investigations of other binary mixtures in order to get 
readings for different noise to side-to-side movement ra- 
tios. This is likely to prove a difficult task and at present 
all we can conclude is that differential pressure segrega- 
tion may play a role in the experiment in addition to the 
differential driving discussed above. 

To demonstrate that different temperatures alone can 
cause segregation we have studied two species with stan- 
dard p parameters but all ra,p = 1- The temperature 
difference is produced by cxp ^ 40 x tr and also setting all 
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FIG. 14: Segregation of cs for standard parameters but with 
Tea; — > 0.05 X ffci: and (Jpx —> 2 X Gpx and no driving of the 
tray base. Thus we get segregation but without coagulation 
at the top and bottom walls. The left and right walls are 
oscillating and prevent coagulation there. Note that the area 
of the system is 1 /4 that of the experimental system, ps are 
coloured black. 




FIG. 15: Coagulation of cs after 96 seconds due only to tem- 
perature difference Tp=1600*Tc. The system is still evolving, 
the groups of cs in the centre will eventually attach to the sides 
of the system. Note that the area of the system is 1/4 that of 
the experimental system, ps are coloured black. 



7 — > 100 X 7 so that the time constants r = m/7 are suffi- 
ciently small that particles remember the temperature of 
their heat baths rather than only the temperature of their 
previous collision partners. This imposed temperature 
difference causes coagulation of the lower temperature 
particles as shown in figure Whilst these parameters 
are very different from the standard ones, they clearly 
show another segregation mechanism in a nonequilibrium 
system and that the heuristic arguments regarding pres- 
sure differences out of equilibrium are valid. 



C. Further Mechanisms 

One further possible mechanism is that the noise of the 
particles is correlated among the particles but differently 
for the cs and the ps. This could cause segregation in 
the same way as the periodic driving since it also would 
produce different collective motions for the two species. 
The correlation could be caused by, for example, the ps 



all changing from sticking to the tray base to slipping at 
the same time in the periodic cycle. Whilst this mech- 
anism is at least plausible, the collective motions due 
to correlated noise and due to periodic driving (which 
would have caused the correlated noise in the first place) 
would not be clearly distinguishable. At the level of this 
phenomenological model which breaks the driving into a 
periodic and a noise component, any such mechanism is, 
therefore, not meaningful. 

It is possible that there are other segregation mech- 
anisms not discussed here, however we believe that we 
have considered the ones most likely to be relevant in 
the experiment. 



V. CONCLUSION 

In this paper we have introduced and numerically stud- 
ied a (relatively) simple phenomenological model of a 
recently reported granular segregation experiment. We 
have measured the same quantities as measured experi- 
mentally and shown that our model reproduces most of 
the features of the experiment, the most important be- 
ing a transition from a mixed to a segregated state as 
the compacity is increased. This behaviour is not a pri- 
ori built into the model - it emerges from the simple rules 
governing the motion of the particles. This is significant 
as it shows that we have a set of basic features neces- 
sary for an explanation of the experimentally observed 
behaviour. 

We then used our model to investigate and identify 
segregation mechanisms and elucidate the experimental 
behaviour. We showed that the transition from mixed to 
segregated state in the model is caused by competition 
between the different driving felt by the cs and ps, which 
acts to cause segregation, and the noise, which acts to 
prevent segregation. We are led to conclude that this is 
also the main mechanism present in the experiment. We 
have also considered and demonstrated segregation due 
to different pressures and shown that it is possible that 
this might play a role in the experiment. The differen- 
tial driving segregation mechanism is applicable to many 
binary driven systems |llLll3l|. 

This work goes some way to explaining the intrigu- 
ing experimental results of 0, 0, Experiments with 
other particle types in conjunction with more accurate 
experimental measurements, particularly with regard to 
particle positions and velocities, should allow more ac- 
curate comparisons with our model and also refinements 
of the model. In particular more accurate values for the 
parameters used. 

Now that we have shown our model to be relevant to 
the experimental system, it is possible to use it further 
to investigate the 'granular statistical mechanics' of this 
type of system. In particular, it may be of use in devel- 
oping and testing theories for agitated granular mixtures 
before attempting the more difficult task of accurately 
comparing with experiments. 
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